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Abstract. We consider the dissipative heat flow and conservative Gross-Pitaevskii dynamics 
associated with the Ginzburg-Landau energy 



\v g u\ 2 , (i-H 2 ) 2 

M 2 4 £ 2 



posed on a Riemannian 2-manifold M endowed with a metric g. In the e — > limit, we show 
the vortices of the solutions to these two problems evolve according to the gradient flow and 
Hamiltonian point- vortex flow respectively, associated with the renormalized energy on A4. 
For the heat flow, we then specialize to the case where M. = S 2 and study the limiting system 
of ODE's and establish an annihilation result. Finally, for the Ginzburg-Landau heat flow on 
S* 2 , we derive some weighted energy identities. 



I. Introduction 

There is a rich and well-developed theory to describe the motion of vortices arising in both 

the heat flow and Schrodinger dynamics associated with the Ginzburg-Landau energy. The 

setting for much of this work has been on a bounded planar domain or all of M 2 , though by 

now there have also been several efforts to understand such flows in M™ for n > 2. In this 

l 
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article, however, we pose these problems on a compact, simply-connected 2-manifold. The 
motivation is to see how the geometry, and in particular the curvature of the underlying surface 
may impact on the motion laws and stability of collections of vortices in these dissipative and 
dispersive settings. 

To this end, we let Ai be a smooth, simply-connected compact 2-manifold without boundary 
equipped with a metric g and then define the Ginzburg-Landau energy E e on Ai for u : Ai — > 

Cby 



|V„d 2 fl-lw 



2\2 



;i.D Et{u) = j^ + ^_^ dVg . 

This is a simplified version of the full energy which includes magnetic effects, originally intro- 
duced by Ginzburg and Landau [H] as a model to describe superconductivity. To E e one can 
then naturally associate the Ginzburg-Landau heat flow 



1.2) 



u\ - A M u £ = ^(l - \u £ \ 2 ) in Ai x R + 
u £ = u £ on Ai x {t = 0}, 

which most efficiently dissipates energy, and the conservative Ginzburg-Landau-Schrodinger 
flow, or Gross-Pitaevskii dynamics 

iu% - A M u £ = - |w e | 2 ) in M x 1+ 



■ B -- B onMx{t = 0}. 



;i.3) 

The latter arises in studies of superfluidity, Bose-Einstein condensation and nonlinear optics. 
Here A^vj denotes the Laplace-Beltrami operator and the positive parameter e corresponds in 
the full Ginzburg-Landau energy to the reciprocal of the so-called Ginzburg-Landau parame- 
ter. In much-though not all-of our analysis, we will be looking in the asymptotic regime where 
e 1. In this regime, it is by now understood ([151 EH I2H [25]) that energetically reasonable 
sequences, bounded in energy on the order of ln(l/e), can possess at most a finite number 
of vortices, that is, zeros carrying non-zero degree. Hence, the analysis of these two infinite 
dimensional flows can be effectively carried out by tracking the motion of a finite number of 
points. 

In the plane, this program consists of showing that in the limit as e — > 0, the role of the 
energy E e in dictating the dynamics is effectively replaced by a so-called renormalized energy 
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W(a, d) dependent on the finite number of vortex locations a = (ai, a 2 , . . . , a n ) and their 
associated degrees d = (di, d 2 , • • • , d n ). The pivotal role of W was first revealed to great effect 
in the stationary planar setting with Dirichlet boundary conditions in [3]. The asymptotic 
motion law for planar vortices of the heat flow (ll.2p . namely 

(1.4) j t a i = -V ai W(z,d) fori = l,2,...,n, 



was first derived in [16] and [19] under a well-prepared initial data assumption up to the first 
time of vortex collisions, and then was extended more recently in the series of papers [1], [5], 
[6]. 

The corresponding system of ODE's governing the asymptotic behavior of Gross-Pitaevskii 
vortices in the plane, 

(1.5) ^a< = ViW(a,d) fori = l,2,...,n, 

was established in [TU [121 [20] and later refined in [IS] . 

A primary goal of the present article is to establish the analogs of (11. 4p and ( 11. 5p in the 
manifold setting. For this, we will appeal to a result of [1], where the author identifies the 
renormalized energy on a Riemannian 2-manifold. Assuming that the manifold is simply- 
connected, compact and without boundary, one can apply the Uniformization Theorem to 
assert the existence of a conformal map h : Ai — > M. 2 |J{oo}, so that the metric g is given by 

(1.6) e 2f (dx\ + dx 2 2 ), 

for some smooth function /. Thus one may identify a vortex a\ G Ai with a point 6, = h(ai) G 
M 2 |J{oo}. Writing b = (pi, b 2 , ■ b n ) with associated degrees d = (di, d 2 , d n ), a result in 
[lj identifies the renormalized energy as 

n 

(1.7) W(b, d):=nJ2 dlfih) - n ^ d t d 3 In \h -b 5 \. 

In Section 3 we derive the asymptotic motion law analogous to (ll.4p for the heat flow (11.21) 
on a simply-connected, compact manifold, valid up to the first time of vortex collisions, under 
an assumption of well-prepared initial data. We follow the basic scheme laid out in [TBI [TH] • 
As in the planar case, one expects the motion to be logarithmically slow (cf. [HI [191 123]), so 
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one must first re-scale in time. Then the key steps are to establish an identity which allows 
one to localize the rate of change of energy about a vortex (Proposition 13.11) and to derive a 
PDE for the energy density itself, Proposition 13.21 In the present setting this entails a new 
term involving the Gaussian curvature of Ai. The main result of this section is the content 
of Theorem 13 .51 

In Section 4 we turn to the vortex law in the dispersive setting on a manifold and generally 
follow the approach of [TTJ H2] . This involves calculating a formula for the evolution of the 
Jacobian of a solution of (11.31) . Proposition 14.11 We then assume that the initial data u e Q is 
almost energy minimizing, cf. (I4.10p . and after a series of results relating the e-limit of the 
solution to the canonical harmonic map on .M, we arrive at our main result, Theorem I4.4[ 
yielding the motion law 

(1.8) d~h = ~(Vj) bi W(b,d) forz = l,2,...,n. 

The evolutions (11.51) and (I1.8P are known as the point-vortex problem, for the plane and a 
manifold, respectively. Alternatively, one can obtain (II .51) from the Euler equations in the 
singular limit where vorticity is concentrated at points. This Hamiltonian system has been 
studied extensively, and we will not attempt to list all references here, but an excellent source 
is [21], where in particular one can find a description of the state of the art with regard to (I1.8P 
on S 2 . Regarding the connection between solutions of (11.31) and (11.81) on S 2 we also mention 
the recent result [13] relating rotating solutions of Gross-Pitaevskii to relative equilibria of 
the point-vortex problem. 

In the final two sections of this paper we return to the Ginzburg-landau heat flow and 
its asymptotic limit. In light of the dissipative nature of this evolution, one expects that 
generically, vortices will tend to annihilate each other. This is particularly the expectation on 
a closed manifold since the total degree of all vortices must be zero, and for example, there 
is no Dirichlet condition as in [3] to force the presence of vortices, nor are we considering any 
applied magnetic field as in P, [10]. Vortex annihilation results in the plane for (11. 21) were first 
established in [2] for any finite e and later the previously mentioned investigations [1], [5], [6] 
carried this out on bounded planar domains in the regime £ C 1 under very mild assumptions 
on the initial data. In [8], the first author addresses the question of whether there can exist 
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stable vortex configurations in the sense of second variation for the energy E e on a closed 
manifold without boundary when e is small and she also presents an annihilation result valid 
for any e for the flow ( 11. 2ft augmented with a Dirichlet condition on a manifold with boundary. 

In Section 5 we consider the special case Ai = S 2 and analyze the corresponding system of 
ODE's 

— ai = -(y g ) a W(&, d) for % = 1, 2, . . . , n. 

We offer a natural definition of how to extend this system past the time of a collision of 
two or more vortices and then working with this definition, we establish in Theorem 15.71 a 
sufficient condition for annihilation of all vortices in finite time involving an assumption of 
initial clustering of all vortices. We also estimate the time it takes for annihilation. 

In Section 6, we remain in the setting of S 2 but turn to the PDE (II. 2p for fixed e, not 
necessarily small. Here we derive in Proposition 16.11 certain weighted energy identities. These 
in particular provide evidence that a similar annihilation result for the heat flow should hold 
under a clustering assumption analogous to the one from Section 5. From these identities 
it follows immediately that any critical point of Ginzburg-Landau on the two-sphere should 
satisfy moment identities suggesting a balanced placement of vortices, cf. Corollary 16.21 

We begin with a section introducing notation and then proceed as outlined above. 

Acknowledgment. The research of both authors was generously supported by NSF grant 
DMS-1101290 and a Simons Foundation Collaboration Grant. 

2. Notation 

Let A4 be a 2- manifold equipped with metric g and let T(Ai) be the tangent bundle of 
M. For X,Y G T{M), f : M ->■ R, u : M ->• R 2 , and v : M -»■ M 2 , we write u = (ui,u 2 ), 
v = (vi,v 2 ) and define the following notation. 

dv g = yj detg\dx 1 A dx 2 - 

V g f: the gradient of / on Ai. 
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Vg /: the skew gradient of / on M. 
Note that V g f and Vjf G T(M). When M = M 2 , we have Vjf = JV P /, where 

'■CO- 

We will also write 

V 3 m := (V ff ui, V 9 w 2 ), Vju := (V^-Ui, V^w 2 ), 
V 9X F: the covariant derivative of F in the direction of X, 
(•, -) g : the inner product of tangent vectors on Ai, 

and 

(Hess f(X),Y) g := {V 9X V g f,Y) g . 

When the inner product involves the gradient of vector- valued functions, we define the notation 
(•, -) g in two cases: 

(V p «, V 9 /) 9 = (V ff /, V 9 u) fl := ((V g f,V gUl ) g , (V g f,V g u 2 } g ) G K 2 , 

(V p «, V g v) g := (V fl «i, V p ui) p + (V^ 2 , V p u 2 ) fl G R. 

We will omit the subscript <? for M = M 2 . 
For G R and A = (A 1: A 2 ) G IR 2 , we define 

A ± = (A 2 ,-A 1 ), n(0) = (cos0,sin0), t(0) = -n 1 (</>). 

Then for any x el 2 , we define 0(x — x ) to be the angular polar coordinate of x centered at 
Xq such that when x ^ x , 

n(0(x - x )) = ~r~ — ^t. 

\X — Xq\ 

When xq = 0, we will simply write n for n(9(x)) and t for t(9(x))). 
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3. Motion Law for Ginzburg-Landau Vortices on a 2-Manifold 

Let M. be a smooth, simply connected compact surface without boundary, and g be a metric 
on M.. We consider the initial value problem 



(3.1) 



u\ - A M u £ = £(1 - \u £ \ 2 ) in M x R + 
u £ = u £ on M x {t = 0}. 

Here for convenience we will associate C with M 2 and consider u : Ai x IR + — >■ M 2 . Note that 
( 13. ip is the heat flow of the Ginzburg-Landau energy defined by 



r |v v\ 2 (1 - \v\ 2 ) 2 r 

(3.2) E £ (u) = l-^-L + 1 ^ 1 J dv g := e £ ( M )^ g . 

As mentioned in the introduction, from the Uniformization Theorem, there is a conformal 
map h : M. — > 1R 2 |J{oo}, so that the metric g is given by 

(3.3) e 2f (dxl + dx 2 2 ), 

for some smooth function /. Thus we may identify points in M. with points in M 2 U{°°}- 
For b = (61,62) •••>6 n ) and d = (d\,d2, ...,d n ) such that bi G 1R 2 , di G Z \ {0} for all z, the 
renormalized energy can be written as [1] 

n 

(3.4) W(b,d) := 7r^rf 2 /(6i) -7tJ2 d i d j ln \ b i ~ b j\- 

i=l ijtj 

To study the dynamics of vortices, we rescale the time variable by a factor |lne| and set 
v e (p,t) = u £ (p, |lne|i). Then v £ solves 



(3.5) 



J ±- l vf-A M v e = $(l-\v e \ 2 ) inMx 



on Mx{t = 0}. 

From the heat flow structure, we easily establish that for v £ solving (13.51) . one has 

d F 

(3.6) - j^e e (v £ )dv g <0, 

which leads to the global existence of a solution to (13. 5p for smooth initial data u £ Q . A key 
tool in capturing the motion of vortices is the following: 
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Proposition 3.1. Let v £ be a solution to (I3.5P and let rj : Ai — > ffi. be any smooth function. 
Then we have 

d_ 

It 



j 7]e £ (v £ )dv g =| lne| / [((Hess r])(V g v £ ), V g v £ ) g - A M r]e £ (v £ )} dv g 
Jm jm 



(3.7) 



V\v £ t \ 2 dv g . 



IM JM 

Proof. Taking the scalar product of both sides of (13.51) with r\v\ and integrating we have 
1 

|lne| JM 



(3.8) 



V\vl\ 2 dv g 



F a e , d f V(v £ 

T)V t ■ A M v dv g = -— / 77— - 

M al JM E 



where V{u) 
(3.9) 



(i-H 2 ) 2 



. Integrating by parts for the second term we obtain 



rjv £ ■ A M v £ dv g 



M 



[ v £ ■ (V g rj, S7 g v £ ) g dv g - — f rj 

JM al JM 



,e|2 



— dVg. 



Combining (13 . 8H and ( 13. 9 p gives 

d_ 

dt 



(3.10) 



7]e e (v £ )dv g 



(VgT], VgV £ )gdVg 



I M JM >M 

To calculate the first term of ( 13.101) . we use ( 13.51) again. Then 



V\v £ t\ 2 dv g . 



■(V g r},VgV £ )gdv g 



M 



(3.11) 



1 



■ (VgV, VgV £ ) g - (Vg V , Vp^) g 



dv„ 



|lne| 



M 



V{v £ 



&MV £ ■ (Vg?7, Vglf)g + — T] 



dv„ 



Through integration by parts and the geometric identity ([22J, p. 207) 



we have 



(V 9 V, y 9WgV sV g v £ )g = (V g v £ ,V gVgV V g v £ ) g , 



[ A M V £ ■ (VgV, y g V £ ) g dVg = - [ (VgV £ , Vg(Vg V , V >') g) 

JM JM 



M 



^IVI 2 ' V 9 r/),+ ((Hess ri)(V g v e ), V g v £ ) g 



dv„ 



(3.12) 



M 



\V g v 



E|2 



A M V - ((Hess 7]){V g v £ ), V g v £ ) g 



dv„ 



Therefore by (l3TT0|) . (EOT]) , and (13TT2|) we derive (l3T7j) . 



□ 
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At this point we assume that the initial data is well-prepared in the following sense. We 
assume there are precisely n zeros {a\, af , a £ n } C M. of u% such that 

(3.13) R := - min min dist(a £ , a £ ) > 0, 

3 0<e<l J 

where dist(p, q) refers to geodesic distance between p and q on M. and 

n 

(3.14) di := deg{u £ , dB 9 R (a £ )) G {1, -1} for 1 < i < n, d, t = 0. 

i=l 

We further assume that for some positive constant c independent of e we have the bounds 

(3.15) / e £ (u £ )dv g < nir\ \ne\ + c, 

Jm 

R 

(3.16) sup {e £ {u £ Q {j))) : dist(p,af) > — for all i = 1,2, ...,n} < c, 

£6(0,1] 2 

(3.17) jnf {\u £ (p)\ : dzs*(p, a?) > | for all i = 1,2, ...,n} > ~, 

(3.18) |«ol < 1, e \^g u o\ +e 2 |Hess u%\ < c. 
We also assume that there exist points C M. such that 

(3.19) limaf = for 1 < i < n, 

£->0 

and 

(3.20) in the sense of distributions. 



i=l 



The last condition in fact comes from the stationary results in [3] (Theorem VII. 2, Theorem 
VII. 3). Next, to obtain the regularity result, we show that the energy density e £ defined in 
(13. 2 p satisfies a certain PDE. 

Proposition 3.2. Let u £ be a solution to (13.11) . Then the energy density e £ (u £ ) satisfies 



[e £ (u £ )] t - A M [e £ (u £ )] =1(1 - \u £ \ 2 )\V g u £ \ 2 ~ |Hess u £ \ 2 - \\u £ ■ V > 



(3.21) - 1(1 - |« £ | 2 ) VP - |VXI 2 ^, 



where Km is the Gaussian curvature of Ai. 
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Proof. From the definition 



Ae 2 



we have 



(3.22) [e £ M] i = (V/,VX) s 

and 



1 — lit 



.= 12 



IV ? /£|2 1 _ | 7 .e|2 
Ajm^tO] = ■ ^MU S 

2 e z 

1 — \v £ \ 2 9 
(3.23) __^J_| V ,«-| 2 + -|««.V^| 2 . 



Substituting (13. ip into (I3.22p we derive 
[e £ (u £ )] t =(V g u £ ,V g 



1 l _ Ue|2 



Aa,u £ + ^(1-|m £ | 2 



)g J 2 U £ ■ A M U S 



4a - i«i a )vi 2 



^4 

1 — |-u e | 2 2 
=(V/, V g A^ m £ ) 9 + — -^\V 9 u £ \ 2 ~^\u £ - V g u £ \ 2 

1 — \v e \ 2 1 

(3.24) l f^u £ ■ A M u £ - -(1 - K| 2 ) V| 2 . 

Using the Bochner-Weitzenbock formula (see [22], Chapter 7, Proposition 33): 



IV u £ \ 2 

A M " ' = |Hess u e | 2 + (V g < V 9 A^ m £ ) 9 + |V 5 u £ | 2 ^, 



and combining (I3.23P and (I3.24p . we have the desired identity. □ 

Since M. is smooth and compact, K M is bounded. Hence there exists e_m > such that 
the proof of the main regularity result from [16] can be readily adapted to the setting on a 
manifold if we set < e < Em- 

Lemma 3.3 (Regularity, cf. [H], Theorem 6.1). Let B®(q) be a geodesic ball in Ai with center 
at q G M. and radius r. Let < e < min{l, r, Sm\- Suppose that u £ is a solution to (13. ip in 
B 2r(l) x (0,4r 2 ) such that 



sup{ / e £ (u £ (-, t))dv g :te[0, 4r 2 ]} < k x . 
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Then there is a constant C = C(k\,M) such that 

e £ (u £ (p, t)) < — for (p, t) G B?{q) x [r 2 , 4r 2 ]. 

If in addition 

e £ (u £ (p,0)) <h forpe B 9 2r (q), 

then 



e £ (u £ (p,t)) < £ for (p,t) G B°(q) x [0,4r 2 ]. 



With Lemma 13.31 we may extend the following result stated in [T7] to our setting on a 
manifold: 

Theorem 3.4 (cf. [17], Lemma 5.1, Proposition 5.5). Suppose v £ is a solution to (13. 5p in 
M. x [0,T ]. T/ien there exists a subsequence {e n } and a finite set of points {cij(t)}™ =1 C M. 
such that Oj(0) = ctj, and /or i? defined in (13.131) . we have 



I e e (v £ (-,t))dv g > ir\n c for all 1 < i < n, e G (0, 1). 



(3.25) 
Moreover, as e n — > 0, 



~n — r ^ 



i=l 

in t/ie sense of distributions for t G [0, Ti) ; where 

T x := inf{t > : min{ dis^a^t), = 0}. 

Let bi : [0,Ti] — > K 2 6e defined as bi(t) = h(ai(t)), and let = Xa^i; where 6i{t) = 
6{x — hit)). Then, after perhaps an adjustment by a constant rotational factor e %e ° , we have 

(3.26) v £n ^v* = n(0) 

uniformly on any compact subset of Q := {(x,t) G IR 2 x [0, Ti) : x ^ bi(t) for all i}. Further- 
more, 

(3.27) \Vv £n \ 2 and2e £n (v £n ) |W| 2 in L} oc {n). 
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Since the method of proving Theorem 13.41 used in [T7] is independent of the geometry, we 
can apply the same proof here with only minor changes. Indeed, for a compact manifold M. 
without boundary, the proof of the convergence is even easier since we do not have to deal 
with the Dirichlet boundary condition. This means, in particular, we do not have to adjust 
the phase of the limiting map v* for specified boundary behavior. We now come to the main 
result. 



Theorem 3.5. The functions {bi(t)} given in Theorem \3.4\ are differentiable, and they satisfy 
the system of ODE's: 



(3.28) 



ik = -l{V g ) bi W{h,d) forte^T,) 
bi(0) = Pi for i = 1,2,..., n, 



where Pi = h((Xi), b = (b±, b 2 , b n ), and d = d 2 , d n ) such that each di is given by 
fl3HD. 

Proof. Fix t e [0, Ti). Without loss of generality, we may set i = 1 and assume that b\(t) = 0. 
For any AeR 2 and p > 0, consider a smooth function rj compactly supported in B 2p such 
that rj = (A, x) in B p . Then for < 5 < p small enough, we have 

riihit + 5)) - 77(0) = (h{t + 5)- h(t), A). 

Now if we integrate 03. 7p from t to t + 8 and divide by | In e\ , then Theorem I3.4I implies 

^(t+5)) - V (0) 

= - lim ( / / [((Hess r/)(V 5 v £ ), V g v £ ) g - A M r]e £ (v £ )] dv g dr 
Ut Jm 



1 



lne| 2 J t 



t+s 



(3.29) --f- / / V \vt\ 2 dv g dr } . 



M 



As a consequence of f]3. TI) with rj — 1, the upper bound for E £ {yl) provided by ( 13.151) and 
the lower bound for E £ (v £ (-, t)) coming from (13.251) . the second integral on the right-hand side 
approaches zero. Since i] is linear in a neighborhood of bi and compactly supported away from 
bi for all i ^ 1, the support of Hess rj does not contain {&i(t), b 2 (T), b n {r)} for r e [t, t + 8] 
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for small 5. Hence by ( I3.26P and ( I3.27P we have 



13 



(3.30) 



t+s 



t JM 



((Hess r)){V g v*),V g v*) g -A M r} 



\V n v 



* 121 



dvgdr. 



Note that for any 77 : M. — > M. and v : M. — > M 2 , using the metric ( 11. 6p we calculate 



(3.31) 



IV v\ 2 

((Hess ri)(V g v), V g v) g - A^Hr 1 " = e" 4/ (gi - 2q 2 + q 3 ) 



where 



(3.32) 



q x = ((D 2 ry)(W), Vv) - Ar] 



\Vv\ 2 
2 ' 



(3.33) 



q 2 = {Vv,Vf)-(Vv,Vr}), 



(3.34) 



?3 = |V^| 2 (V/,V?7}. 



Now combining ( 13 . 3 1 [) - (13 . 34 j) we obtain 



M 



((Hess r/)(V/), V g v*) g - A M r] 



\V g v*\ 2 



dv„ 



-2/ 



B 2p \B p 



((D 2 r}){Vv*),Vv*) - At] 



*|2 



dx 



(3.35) 



+ / e" 2/ [|W| 2 (V/, V77) - 2(W, V/) • (W, Vtj)] dx. 

JB 2p \Bp 
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Let n = (rii, riz) be the outward unit normal. Integrating by parts we derive 



-2/ 



B 2p \B p 



((D 2 ??)(Vt;*),Vt;*) - A v 



B 2p \B p 



iwr 


dx 


2 




1 * 12 
IT 
1 u x 2 1 



y a;i v x 2 )'\x\x 2 



dx 



B 2p \B p 



e 2f [Av* ■ {v xi v* xi + Vx 2 v* X2 )] dx 



'B 2p \Bp 



-2/ 



xi J x 2 r lx 2 ){\ v x 1 



* |2 L,* |2 



<J 2 ) + 2(/ xl »7sa + fx 2 Vx 1 )vl 1 ■ v* X2 ] dx 



-2/ 



(^"l - W^KKJ 2 - KJ 2 ) + fe 2 ^i + Vx 1 n 2 )v* 1 ■ v* 2 



ds 



'B 2 p\B p 

(3.36) + / e 

J dBp 



[ e~ 2f [-|W| 2 (V/, Vrj) + 2(W,V/) • (W,Vr/)] 

JBorXBn 



-2/ 



|Wf (Vr/,n) - (W\n) • (W\ Vr/) 



The last equality comes from the fact that v* is harmonic in Bi p \ -Bp. Then by (I3.35f) and 
(I3.36P we have 



M L 



((Hess r?)(V^*), V^*) ff - A mV 



\V g v 



* 12 





}|V.'| 


'dBp 





■|WT(V?7, n) - (Vu*, n) • (W\ Vr/) 
(3.37) :=/i(r)+/ 2 (r). 

Recalling from Theorem 13.41 that v* = n(O), we evaluate the following terms at x = pn: 



|W| 2 n 



4 + ^ <Ed W t) 



n + Ei(x, r) 



(3.38) 



i 9^7 
P P 



n + Ei(x,r) 



(3.39) 



(W, n) • \7v* = —(C(x, r), n)t + E 2 (x, r), 
P 
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where 

(3-40) C (^) : =X>ur% 

1=2' 1 ' 

and Ei, E2 are bounded. Next we expand e~ 2 ^ near x = 0, 

(3.41) e~ 2 f = e~ 2 ^ - 2 e - 2 ^ (V/(0), pn) + 0(p 2 ). 

Since the most singular part, e~ 2 ^°)-^-n, from (13.381) makes no contribution to the integral in 
f l3~37D . substituting (13T38D . ( 13391) . and (131TD into Ii(r) and J 2 (r) we obtain 

(3.42) h{r) = e- 2 ^ [*" [rf 1 (C(x, r), t) - (V/(0), n>] (n, V V )d6 + O r (p), 

Jo 

and 

(3.43) J 2 (r) = -e" 2 ^ / * di(C(a:, r), n)(t, V V )d9 + O r (p), 

Jo 

where T (p) indicates a quantity which is O(p) with the implicit constant depending only on 
t. Observe that for any fixed V e I 2 we have 



/ (V, t)nd6 = ttV^ = - / (V, n)td0, 
Jo Jo 



and 

-27T 



\ (V, n)nd# = ttV. 



Thus 

I^r) + J 2 (r) = -Tre" 2 * * (V/(0) - 2d 1 C- L (0, r), Vr/(0)} + O r (p) 

(3.44) = -7re- 2/ W(V/(0) -2d 1 C" L (0, r), A) +O r (p). 

Combining 1ET5D)) . (13T3TD . and (EH we derive 

= (lim-e- 2 ^- / [V/(0) - 2d 1 C ± (0,r) + T (p)]dr, A) 

(3.45) = (-e- 2 ^(*»[V/(6i(t)) - 2d 1 C- L (6 1 (t),t)], A> +0(p). 
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Note that (ET351) holds for all A G M 2 . Also since 

n 

W{b, d) := 7T f(bi) - 7T d t d 3 In \h -b,\, 

i=l i/j 



we have 



-(V 9 ) bl W(h,d)=e 



-2/(61) 



7T 



V/(6i)-2di53d j 



i=2 



h - bj 



e -2/( 6l ) [v/fa) - 2d 1 C d -(6 1 )] 



Then taking p — > in ( I3.45P we obtain 
(3.46) 



dh 1 



- z (V fl ) 6l W(b,d). 



4. Motion Law for Gross-Pitaevskii Vortices on 2-Manifold 
In this section we consider the initial value problem 



(4.1) 



%n% - A M u £ = 4(1 - \u £ \ 2 ) in M x 



on Mx{t = 0}. 



□ 



for u : M. — > C. Global in time well-posedness in this defocusing setting is provided, for 
example, by the results in [7] 
We define the current 



j(u) := (iu) ■ W g u. 

Then the relationship J(u) = ~V x j(u) between the Jacobian of u and the current on R 2 
motivates the following definition of the signed weak Jacobian J(u) for u G H l (A4) via 

(J(u),v) = I I (j(u),Vjr)) g dv g , for all 77 G C\M). 

Note that J(u) can be viewed as an element of the dual of C 1 . We first establish the following 
identity, which will be used in our study of dynamics of Gross-Pitaevskii vortices. 
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Proposition 4.1. Let u £ be a solution to (14. ip and let rj : Ai M. be any smooth function. 
Then we have 

(4-2) J t {J{u£) ^ = ~L {VsUS ' V ^ V ^M*V 

Proof. By the definition of J(w), we have 

Tt {J{ue) > v) = \S M { Jt j{u% v iri° dv °- 

Applying ( 14. ip we deduce 



j t 3(u £ ) =(iu £ t ) ■ S7 g u £ + (iu £ ) ■ WgU 



[A M u £ + ^{l-\u £ \ 2 )u e )-V g u 



- u £ ■ V g {A M u £ + ^(1 - \u £ \ 2 )u £ " 



Thus 
(4.3) 
where 



and 



j t (Au £ ),r ] )= 1 -(A 1 + A 2 ), 



A x = I {& M u e + \{l-\u e \ 2 )u e )-{V g u e ,V^ri) g dv g 
Jm e 

V(u £ ) 

A M U £ ■ (VgU £ , V^), - (V,-^, Vj V )g 



dv 



A 2 = - I u £ ■ (V g (A M u £ + -1(1 - \u £ \ 2 )u £ ), Vjr])gdv g . 
Jm £ 



Note that divV^r] = 0. Then integrating by parts we obtain 



At = / A M u £ ■ (V g u £ i Vfv)gdv B 
Jm 



(VgU £ ,Vg(VgU £ ,Vj V )gdVg 



M 



M 



(V 9 « £ , V^Vtv), + ^(V g | V g v £ \ 2 , V>) 



g 'II g 



dv n 



(4.4) 



Jm 
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and 

A 2 = f (A M u e + \(1 - |w £ |V) • diy(u £ Vjrj)dv g 

J M e 

= [ (A M u £ + \(1 - \u £ \ 2 )u e ) ■ (V g u e ,Vj V } 9 dv 9 

J M £ 

(4.5) =A,. 

Combining (14. 3p . (14. 4p . and (14. 5p we have the desired equality. □ 

For {m £ }o< £ <i C ^(M) satisfying the almost energy minimizing assumption (14. 6 p given 
below, the following upper bound results are proved in [12] for sequences {J{u £ )} converging 
weakly as measures. A similar result is also established in [2"0] . Proposition 3.3. Since the 
proof is independent of the geometry, the results adapt without change to our setting on a 
manifold: 

Lemma 4.2 (cf. p2], Theorem 1.4.5). Let {v 6 } C H l (M) be any sequence such that 

n 

J(v e ) -»« n^didai 
i=i 

weakly as measures. Here a« G Ai, and di G {±1} with *YTi=\ di = 0. Suppose that there exists 
some 7 > such that as e — > 0, 

(4.6) E £ {v e ) < nn\ lne| + W(b, d) + 7 + o(l), 

where d = (di,d 2 , ...,d n ), b = (bi,b 2 , ...,b n ), and bi = h(ai) (cf. (Il.6p - fll.7p ). Then there exists 
a constant C such that for every p > 0, 

( 4 - 7 ) limsup £ ^ ||^ - j(H)\\ 2 L2{MX[J7=iBap(ai)) < Cj, 

( 4 - 8 ) limsu P^ollv 9 |^ini2 {Mur=ii?pS(ai)) <C 7 . 

Here H = H(sl, d) is the canonical harmonic map, which is unique up to an arbitrary rotation. 
By canonical harmonic map we mean that H(&, d) is a harmonic map into S 1 with singularities 
at points a = (ai, a 2 , a n ) with Oj G M. such that the winding number (degree) of H about 
ai is di with d = (di,d 2 , d n ). 
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Now we make the following assumption on the initial data u e Q . We assume that f)3. 13j) . 
(13.14p . (I3.15p . and (13 . 191) hold. Furthermore, we assume that 



where is given by the assumption (I3.19p . Finally we assume that u e Q is almost energy 
minimizing, i.e. for every p > 0, Uq satisfies 



as e — )■ 0, where (3 = (fa, fa, Pn), and Pi = h(a>i). With these assumptions, the following 
results established in [12] carry over without change to our setting. 

Theorem 4.3 (cf. [12], Theorem 1.4.1). Suppose u e is a solution to (14. ip with the initial data 
Uq satisfying the above assumptions. Then there exists a subsequence {e n } — > and a finite set 
of points {ai(t)}^ =l C Ai such that aj(0) = functions {aj(t)}™ =1 are Lipschitz continuous, 
and 



n 



(4.9) 




i=l 



(4.10) 




n 





i=l 



fort G [0, Ti), where 



(4.12) 



T\ := inf{t > : min{ dist(cii(t), dj(t))} = 0}. 



Moreover, 



(4.13) 



u £n \ 1 zn L 2 ([0,T!];L 2 (^)) 



and 



71 



(4.14) 




i=l 



where H = H(a.(t),d) is the canonical harmonic map. 



We can now establish our main result of this section. 
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Theorem 4.4. Let functions {aj(t)}™ =1 be given as in Theorem \4-3[ Then for bi(t) = h(a,i(t)), 
the collection {&j(t)}™ =1 C IR 2 satisfies the system of ODE's: 



(4.15) 



<kibi = -${Vi) ii W{b,d) / O rte(0,T a ) 
6i(0) = 



where 0i = h(ai), b = 6 2 , ■ •-, &n)> d — (^i, <^2, ■ •-, <^n)- 
Proof. -Step 1: 

Let the functions {bi(t)} solve the system of ODE's (14. 15j) let T\ be the time of existence 
so that 

fx := inf{t > : min{|&<(t) - ^(t))} = 0}. 

Define the function £ by 

n 

C(t) = J2\ d Mt)-dMt)\ 

8=1 

for t G [0,2;) where T» := min-pi,^} and 7\ is given by fH7l2]) . Fix a > 0. Since C(0) = 0, 
there exists < t a < T* such that ((t) < cr for t G [0,t CT ]. Since the functions 6j(t) are 
Lipschitz continuous, they are differentiable for almost every t. Hence, for a.e. t G [0, t a ], we 
have 

— < \d — d ^1 
dt~ ^ 1 dt i dt 1 

^ E 1^ + -W^(b, d)| + - J2 l(v^)^(b, d) - (vjy b w(b, d)| 

i=l i=l 

(4.16) < £ \d~ + -(Vj) bt W(b, d)| + C(. 



i=l 



Here we have applied Taylor's Theorem to W, valid for a sufficiently small. 

Now, fix t G [0, to be any point of differentiability of 6j(t). Without loss of generality, 
we may set % = 1 and assume that &i(t) = 0. Pick A G IR 2 such that 

(4-17) \d^ + i(V^) 6l W(b,d)| = (d^ + ±(Vj) bl W(b,d), A). 
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Fix p > and let 77 be a smooth function compactly supported in B 2p such that r\ = (A, x) 
in £> p . Then applying Proposition 14.11 and (14. lip we deduce 

(d! ^ (t) , A) = (lim d -f [h (t + 5) - h (f )] , A) 



1 r f 

(4.18) =- lim lim— / / (V q u e , V„ w e Vtrj) Q dv q dr. 

-Step 2: 



Let if be given in Theorem 14.31 Projecting if onto the plane, we may write H = u* := n(@), 
where 6 = Y^i=i and 9i(t) = 9(x — We will show that 

(4-19) ((\7j) bl W(h,d),A)= I {V g u\V 9Vgu ^r,) g dv g . 

JM 

Note that for any smooth r] : Ai — > R and it : .M — )■ M 2 in if 1 , using the metric (13. 3p we 
calculate 

(4.20) - (V,u, V 5Vs „V» 9 = e" 4 % + 2q 2 - g 3 ) 



where 

(4.21) qi = (Vu, (£> 2 t?)JV«) with J = 

(4.22) fcHV^V/MV^V-S?), 
and 

(4.23) g 3 = |V M | 2 (V/,V ± r/). 
Taking < r < p, and combining (I4.20l) - (l4.23p we derive 

- / (v/,v 5v v 9 ^)A 



1 
-1 



e- 2f (Vu*, (D 2 r/)JVu*)dx 



(4.24) 



/ e~ 2/ [2(Vit*, V/) ■ (W, V^) - | Wf(V/, V^r/)] dar. 

JB 2 r\B r 
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Integrating by parts and using the fact that u* is harmonic in B 2r \ B r we obtain 



I e~ 2f (Vu*, (D 2 ri)Wu*)dx 

J B2r\B r 



B2r\B r 



a 



U xi\ )Vx 1 x 2 + ( U Xl ' U x 2 )( r l x i x i VX2X2)} dx 



2/ \(\ u * |2 j_ |„* |2 



B2r\B r 



XI I 



+ WZaWiVx* - (KJ 2 + K 2 l 2 k^i] dx 



+ 2 / *- 2 / n«* i 2 



e ' [|< a | fx 2 Vxx ~ |u*J /x^x 2 + «! • u* X2 )(f xl rj xl - f X2 r) X2 )\ dx 



+ / e 2/ [l<| 2 % 2 n i - + « x ■ u* X2 )(r] X2 n2 - r) xi n x )\ ds. 

J a B r 



Next, integrating by parts again for the first integral, we have 



e~ 2f (Vu*, (B 2 r])Wu*)da 



B2r\B r 



B 2r \B r 



dB r 



e 2/ [(KJ 2 - K 1 \ 2 )(fx2Vx 1 + fxiVxa) + 2« • u* X2 )(f Xl r] Xl - f X2 rix 2 )] dx 



;(KJ 2 - KJ 2 )fein 2 + rfe^i) + (u* xi ■ u* X2 )(r] xl ni - r] X2 n 2 ) 
•■ ! e~ 2f [|Vw*| 2 (V/, V^) - 2(Vw*, V f) • (Vu*, V^)] dx 

J Bo r \B r 



ds 



(4.25) 



dB r 



-2/ 



|Vu*nV^,n) - (Vu*, V^) • (Vu*,n) 



ds. 



Thus from (I4.24p and (|4.25|) . we have 



(4.26) 



-2/ 



|Vu*nV x r/,n) - (Vu*, V x r?) ■ (Vu*, n) 



(is. 
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Note that V -1 // = A 1 - at x = and that (A- 1 )- 1 = —A. From the type of argument used for 
calculating ( 13.371) . we obtain 



/ (v/,v Sv ,v^)A 

Jm 



IM 

= (e- 2 ^ t ^f(h(t))-2d 1 C ± (h(t),t)] ± ,A)+0(r) 
(4.27) =((Vj) bl W(b,d),A) + 0(r). 

Since r can be taken arbitrarily small, we have proved (14. 19ft . Using it we derive 

-t+8 



(-(Vj) bl W(b(t),d),A) = limj- [ ((Vj) bl W(b(r),d),A)dr 
(4-28) = hmj-J^Jj\/ 9 u*, V 9VgU *Vjv) 9 dv 9 . 



-Step 3: 

From ( HTTP . KW and (14281) . 



dbi 1 
dt 7i 

-t+S 



(4.29) = - limlim j- J jj V S u ^ V ^ V» s - (V 9 u*, V 9v ^ V^^dr. 
Applying ( 14^201) -( 14T231 ) again we find 

= / e" 4/ [(Vu £ , (£ 2 t/)JV?/) - (Vu*, (D 2 r])Wu*)}dx 
Jb 2p \b p 

+ 2 [ e" 4/ [(Vu £ , V/) • (Vu £ , V ± r/) - (Vu*, V/) • (Vu*, V ± r/)]rfx 

(4.30) -/ e- 4/ (|Vu £ | 2 - \Vu*\ 2 )(Vf,Vr])dx. 

Jb 2p \b p 

Note that in (I4.30p . each term in the integral is of the form Fu e x u%, where F is a smooth 
function and k, I G {1, 2}. To control them, we first calculate for k = 1, 2, 



£ j fc (u e ) iu e u e 

Uxk ~ \u £ \ W\ + ' M ' Xfe Wv 
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where we write j{u £ ) = (j 1 {u £ ) , j 2 (u £ )) . Thus for k, I G {1,2}, 
j k (u e )j l (u £ ) *w, 



\u 



e|2 



-r(0 



(4.31) 



+ /(«*! 



J " ( " E) -/(«■•: 



|u e | 



Since is uniformly bounded in L 2 (Ai \ UILi ^p( a i) x + ^])> U P to a subsequence it 
converges to some weak limit j in L 2 . Then from (14. 13j) and (I4.14p . we see that 

3 =j(u*), 



i.e. 



(4.32) 



j(u £ 



j(u*) in L 2 (M\\jB^ ai )x[t,t + 5]). 



i=l 



Now the almost minimizing energy assumption (14.101) and the conservation of energy E e along 
the flow (14. ip imply 

< mr\lne\ + W((3,d) + o(l) 
= nn\\ne\ + W{b{t), d) + o(l) 
(4.33) < mr | lne| + W(b(t), d) + CC(t) + o(l). 

Here we also used the fact that the renormalized energy W is conserved by the flow b(£). 
Then by (14. lip and Lemma [4.21 we have 



(4.34) 
and 
(4.35) 



lim sup 



limsup £ ^ ||V g |-u £ ||| L2(A4VU n =i ^ (ai)) < C((t). 



Using ( 14.32(1 and passing limit e — > in ( 14. 29(1 we obtain 

\di^ + l(yi)mW(b,d)\<CC{t). 



GINZBURG-LANDAU AND GROSS-PITAEVSKII VORTICES 25 

Since such an estimate holds for each vortex path we combining this with ( I4.16P to see that 
for t € [0,£ CT ], one has 

with ((0) = 0. Hence, we have ( = 0, which implies (I4.15p . □ 

5. Gradient Flow of Point Vortices on S 2 

In this section we will discuss the limiting vortex motion on S 2 C M. 3 given by the system of 
ODE's (ET2g|) . Let {P,}^ be 2n vortices on S 2 with degrees {di} 2 ^ such that Y%Li d i = °- 
Here we restrict attention to vortices of initial degrees ±1. Note that for S 2 , we may write 
the renormalized energy W as 



(5.1) W = -TT^^lnlPi-P 



Defining T\ := inf{t > : minj^{|Pj — Pj|} = 0}, the vortices {Pj} 2 =i satisfy the gradient 
flow: 

dP 

(5.2) -^ i = (V sa ) Pi diX^d i ln|P < -P i | R 3 forte(0,T ). 

Now let Pi G R 2 be the image of Pj via stereographic projection such that coordinates of Pj 
in R 3 are given by 

(5-3) P < = (T^>fe4). 

1 + rf 1 + rf 

where rf = \pi\ 2 . In this case the conformal factor e 2 ^ x ^ = ,,.^, 2 \2 , i-e. 



f(x) = ln(- 



(l+|x| 2 ) 2 

2 , 



12'" 



1 + \x\ 

Thus (I3.28P can be written as the following system of ODE's on the plane: 
dpi _ . (1 + r 2 ) 2 / p, ^ ^ - pj 



(5.4) ^P = p\ = + 

We first establish a useful identity for the quantity Y^l=i 
Proposition 5.1. Let V = Y^Zi P*> then V = V in (0, Ti). 



2 
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Proof. Taking the derivative using (I5.3P we find 



2» 



(5.5) 



Vo = (£ 



i=l 



1 + r, 



2 P* 



4pi ■ Pi 



(1 + r 



2)i 



E 



4pi • pi 



r(i+^ 2 ) 



2N2^ 



Then using ( 15.4ft we obtain 
2 4 Pi • Pi 



2N2 



P) 



P» + + r}) J2 dj | P * Pj , 2 - Pi [ T~T~2 + 2 ^ E ^ 
^ - Pil 2 \1 + ^ 2 ^ |p.)-p 



2r? 



r i - Pi • Pi 



■ |2 



P ' 1 + rf IP.-P.I 2 J ^IP^P,! 2 



(5.6) 



i^i 



Pi ~ Pj 
3l Pi~Pj\ 2 



Since rf — 2pi • Pj = |pi — Pj| 2 — t 2 and di Ylj^i dj = — 1, we have 



2 . 4pi ■ p^ 

i , 2 P* /i i 2N9Pi 

1 + rf (1 + rf) 2 

2r 2 
1 + rf 



E 



iPi -P 



Pi ~ Pj 
Pi - Pi 



(5.7) 

Similarly, 

(5.8) 

where 



2pi 



1 + rf 



A,, 



4pi • pi rf - 1 



1 + rj 



(1 + r. 



2N2 



5j = C?iC?j- 



jVi 



)>>2 

i i 

Pi - Pil ; 



Note that Y%L X A * = and X]-=i B { = 0. Combining (E2J), (EZD and (EH]), we derive 



^2)i 



V = V . 



□ 
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Corollary 5.2. // V (0) ^ 0, then 

dv ° 0, 



dt |V | 

i.e. the projection of V onto S* 2 zs fixed for t e [0, Ti). We may mew i/iis as conservation 
of the center of mass. Furthermore, either V (t) = or else Ti < oo. Lastly, \ = is a 
necessary condition for an equilibrium solution to (15. 2ft . 

Proof. By Proposition 15.11 we have V (t) = V (0)e* and since |V (t)| < 2n, all of the 
conclusions follow immediately. □ 

Proposition 5.3. Let V be as in Proposition \5.1[ Then Yli>j ~ F*j | 2 decays in [0,7i). 
Furthermore, the decay is strict z/Vo(0) 7^ 0. 

Proof. Since |Pj| = 1 for all i, we have 

l>j l>j 

= -|(v„.v„) 

= -2V ■ V„. 

From Proposition 15.11 we then have 

l^lP.-P/^^lVoi^^lVo^lV* 

i>j 



□ 



If Tx < 00, we know that = Pj t at t = T± for some z* 7^ j*. Next we will prove that any 
collision must involve at least two vortices with different signs. In general, we have 

Proposition 5.4. Assume that at the first collision time t = T\ < 00, the total degree of the 
colliding vortices is not zero. Denote the absolute value of this total by I. Then the collision 
cannot involve only vortices having degree of one sign. Furthermore, if the collision involves 
k vortices of one sign and k + I vortices of the opposite sign with I > 1, then the following 
inequality gives an upper bound on I: 

C k 2 + C\ M - k{k + /)<!, 
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where C„ 



ki 



n\(k—n)\ ' 



Proof. Without loss of generality, we assume a collision occurs at the south pole, which in our 
coordinates on IR 2 corresponds to the origin. We will consider the case where the number of 
colliding vortices with degree 1 exceeds the number with degree —1. The other case is handled 
similarly. 

Thus we let k denote the number of vortices with degree —1 and let k + l denote the number 
of vortices with degree 1. We also denote by / the set of indices i G {1, 2, 2n} corresponding 
to vortices involved in the collision at origin. Using (15. 4p we calculate 

^ r 2 _ _^ — ft i ^2 { r i , A A r i ~ P * ' Pj 



2 P , l -p l = (l + r^[ T f-, + d l J2d 

By the definition of / 



lim pi(t) ^ for % ^ /. 

t— >T\ 



Thus for z G / and < T\ — t <C 1, 



Therefore, 



E^ 2 = E E ^ L_p.p +°(i)= E 



Note that in / there are A; +1 vortices with degree 1 and k vortices with degree —1. We easily 
check that 

E didj = 2[C k + C k+l - k(k + I)} > 2, 
where the last inequality follows from contradiction hypothesis: 
(5.9) C k + C k+l - k(k + l)>l, 

and we use the convention that C° = C\ — 0. Note that (15. 9ft is satisfied if k = since in this 
case I > 2. But this implies that there exists z* G / such that ^r 4 2 > ^7 for all t sufficiently 
close to Ti, which contradicts the assumption that P j — > as t — > T\ for all % G I. □ 
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When two or more vortices collide at some time T , the system (I5.4p is no longer well- 
defined. However, we wish to consider a natural continuation of the flow after this collision 
time by restarting the flow with an appropriate removal or recombination of the colliding 
vortices. This motivates the definition below. We first make the assumption that if for some 
P* G § 2 and some / C {i 6 N : 1 < i < 2n}, vortices {Pj}j 6 / collide at P* when t = T*, we 
must have 



Remark 5.5. When there are 2 vortices, it is clear that assumption (I5.10P holds. For 4 
vortices, Proposition 15.41 indicates that collisions cannot involve only vortices with degree 1 
or only —1. Thus assumption (I5.10p is also satisfied in this case. 

Definition 5.6. We assume that if vortices {Pi}j £ / collide at time t = T* ; then they annihilate 
each other ifYliei^i = 0. In this case, we restart (15 .4p for the remaining vortices with initial 
time t = T* and {Pj}j g / dropped. If^2 ieI di = l(resp. — 1), then we assume the colliding 
vortices {Pj}j 6 / combine to form a single vortex at P* with degree d* = l(resp. — 1). Then 
we restart (15 .4p at time t = T* for the remaining vortices with {P,,}j e / dropped and P* added. 

Using the definition above to extend the flow past collisions we may prove the following 
theorem that gives a clustering condition on initial data at t — that implies all vortices will 
eventually be annihilated. Here, for convenience, we assume the degrees {dj}^ satisfy 



Theorem 5.7. For fixed < s < 1, let A s = {(x,y,z) E E> 2 : z < —yl — s 2 }. Assume that 
vortices {Pi} 2 ^ satisfy (15.21) with Pj G A s att = for all i. Furthermore, assume that (15.101) 
holds at each collision time with the flow defined past collisions via Definition \5.b\ Then there 



(5.10) 




(5.11) 




are no vortices after a finite time if \/l — s 2 > 



n-l 



n 



Proof. Let 



V :=Pi + P 2 + ... + P 2n . 
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Since Pj(0) G A for all i and s < 1, |V (0)| > 2ny/l — s 2 > . Then from Proposition 15. 1[ 
since Vo(t) = V (0)e*, there exists the first collision time < T\ < oo. Assume that there 
exists a sequence of m collision times < T\ < T 2 < ... < T m < oo with m > 1. Note that 
since each collision results in the annihilation of at least two vortices, m < n. 

Since (I5.10P holds at each collision time and the total degree of vortices is zero for all t > 0, 
there must always remain an even number of vortices in the restarted system of ODE's (15 Ah 
after the collision. Similarly, in light of Definition 15.61 and assumption ( I5.10p . we note that 
after the collision time Tj., an even number of vortices will be removed. For k = 1,2, ...,m, 
we define the total number of the removed vortices from t = to t = 7\. by 2j^- Recall that 
from (15. lip , the sign of the degree of a vortex Pj depends on the parity of the index i. Thus 
for each 1 < k < m we may relabel the indices of vortices such that at t = 

while 

Pi ^ for i ^ i', i, i' > 2j k , 
where = jo < ji < J2 < ••• < jm < n. Furthermore, for t < Tk 

P2j fc -1 7^ ^2j k - 

If jm — n i then at t = T m all vortices have been annihilated. Therefore we will proceed by 
contradiction and assume that j m < n — 1. Now for each 1 < k < m, set 

Vfc := P2j ft +1 + P2j fc +2 + ••■ + P2n, 

which is the sum of surviving vortices for t > T^. We view each as being defined on 
[Tfc,T/ c+1 ) with Vfc(Tfc) being a sum of vortices after appropriate removal of colliding vortices 
as in Definition 15.61 Then using the definition of V and V^, we have 

2jl 

VoCZl) = U 1 + V 1 (T 1 ), where V 1 := ^ P i( T i)> 

i=i 

Vi(T a ) = U 2 + V 2 (T 2 ), where U 2 := ^ P f (T 2 ), 

i=2ji + l 
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and so on, 

V m _i(T m ) = U m + V m (T m ), where U m := ^ p i( T ™)- 

i=2j m _i+l 

Here for each k we evaluated Vfc_i(Tfc) by taking the limit t — > . Note that after each 
collision time T k , the restarted system of ODE's (15.41) for remaining vortices {Pi}f=2j k +i nas 
the same structure as the original one i.e. Proposition 15.11 and Proposition 15.31 hold for the 



restarted system. Hence from Proposition 15.11 we have 

(5.12) jY k ^ = Y k ^ for t E (T k _ u T k ) 
for 1 < k < m. Applying (I5.12p . we derive 

\T (rp \JT m -T m -\ T T 

v m— ly 1 - m— 1 

— [Y m -2(T m -i) — U m _i]e Tm Tm - 1 — JJ m 

\r (T \„T m —T m -2 tt „T m -T m _i tt 

— v m _2( i i m -2je — u m _ie — u m . 
Continuing this process finally we obtain 

m 

V m (T m ) = VoCTOe" 1 - ^U fc e Mfc 

fc=i 

m 

(5.13) =e T ™[V (0)-^U fc e- Tfc ]. 

fc=i 



Since v 7 ! 3 ^ 1 > — 

| ^ U k e~ Tk | < 2j m < 2(n - 1) < 2nVl - s 2 < |V 



m 



fc=l 



Thus |V m (T m )| ^ 0. But then from Proposition 15.31 there is at least one more collision after 
t = T m which contradicts the assumption that there are only m collision times. □ 



Corollary 5.8. Under the assumption of Theorem \5.7\ all vortices have been annihilated by 
the time T = In -, where 



k := n 



s/l - s 2 - (n - 1) > 0. 



32 K. CHEN AND P. STERNBERG 

Proof. Suppose T m is the last collision time, then 1 < j m _i < n — 1. From ( I5.12p and (I5.13P 
we have 

(5.14) V m _x(T m ) = Vm-i{Tm-i)e Tm - Tm -\ 
and 

m— 1 

(5.15) V m ^(T m ^) = e Tm - 1 [V Q (0) - ^U fe e" Tfc ]. 

fe=i 



By the assumption of Theorem 15.71 |Vo(0)| > 2nVl — s 2 . On the other hand, 

m— 1 m— 1 

|$^U fc e- Tfc | < ^ |U fc | < 2j m _ 1 . 

k=l k=l 

Thus 

m— 1 



it=i 



> e^" 1 (2nVl - s 2 - 2j m _x) 
(5.16) =2e T - 1 («: + r2- l-j m _i). 

Substituting ( 15. 16ft into ( 15.14ft we obtain 

\V m -i(T m )\ > 2e Tm (K + n-l-j m ^). 



Noting that |V m _i| < 2(n — j m -\) and k < 1 we have 

K + n — l—j m _iJ ~ K 



(5.17) T m < In ( ■ 71 { m ~\ ) < In-. 



□ 



Remark 5.9. From Remark 15. 5[ Theorem 15.71 indicates that all Pj's will vanish after a finite 
time if all P«(0) G A s for any s < 1 when n = 1, and any s < ^ when n = 2. In fact, we 
have an explicit solution of (15.41) when n — 1: pi = (g, 0), p 2 = (— g, 0), where 
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and 

c = if q(0) = 1 
c> if g(0) > 1 
c<0 if g(0) < 1. 

Note that g(0) > 1 (resp. g(0) < 1) implies that Pi(0) and P2(0) are in the upper hemisphere 
(resp. lower hemisphere). Both cases lead to collision as predicted in Theorem 15.71 When 
n — 2, after the first collision time T 1; either all vortices are annihilated or there remain two 
vortices having different signs. In the latter scenario, the two-vortex flow then proceeds as 
described above. 



6. Weighted Energy Identities for Ginzburg-Landau on S 2 

In this section we return to the PDE setting of (13.11) and we set Ai = S 2 . As a gradient 
flow, we of course know that solutions satisfy the standard dissipation rule 



\uA dv, 



.</• 



M 



but our goal in this section is to derive weighted energy dissipation rules that we believe 
should have implications for vortex evolution on the 2-sphere. As a by-product, we will 
derive necessary conditions on equilibrium solutions to (13. ip or (14. ip that echo the symmetry 
requirements on vortex placement implied by Corollary l5.21 We wish to emphasize that, unlike 
the analysis in Sections 3 and 4, these results are not asymptotic in e but hold rather for any 
positive fixed e. 

C we define the weighted energies F\, F 2 and F 3 by 

1 — X\ 



To this end, for «:S 2 xl + 
F 1 {u) := 

F 2 (u) := 

F 3 (u) :-- 



S 2 



S 2 



s 2 



\V g u\ 2 



\V g u\ 2 



+ V(u) 
+ V(u) 



2 

l — x 2 
1 -x 3 



dv. 



dv g and 



dv, 



y 



where V(u) = - ^ and we write the coordinates of any P G 5* 2 as P = (x\, X2,x 3 ). Again, 
since the results to follow hold for any positive e, we suppress the dependence of e in the 
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notation for the weighted energies and in the solution u to ( 13. ip . Note that if p G M. 2 is the 
image of P via the stereographic projection mapping the north pole to oo, then we have 

r 2 - 1 

(6.1) x 3 = where r 2 = Ipl . 

1 + r 2 

Transforming u via stereographic project (and still denoting it by u) we note that for example 
F 3 can be written as 



where dx = dx\ dx2, w(s) = and r 2 = |p| 2 . 

Proposition 6.1. Let u be a solution to (13. ip . Then for i = 1,2 and 3, Fi(u) satisfies 



Proof. For simplicity of notation only in what follows we set e = 1. We will derive the identity 
for F3. The other derivations are identical. Taking the derivative of (16. 2 p we obtain 



(6.2) 




(6.3) 





Integrating by parts we have 





Thus 



d 
dt 




Au + (l-\u\ 2 )u + u t ■ (Vm, Vw(r 2 )) dx 



(6.4) 




Here in the last equality we used the fact that from (13. ip . u satisfies 



(6.5) 



ik = 




Au + {1- \u\ 2 )u in R 2 . 



4 
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Since Vtw(r 2 ) = 2w'(r 2 )x = — jr^pzy. , applying (I6.5P again we have 



/ u t ■ (Vw, Vw(r 2 )) dx 
Jm 2 



„2\2 



-Att + (1 — kt )w 



4 

;i + r 2 ) 2 

4 



Am + (1 - |w| 2 )w 



2k; (r )(Vu, x) (ix 
2 



;i + r 2 ) 2 



(Vw, x) dx. 



Then integrating by parts twice we find 



/ u t ■ (Vu,Vw(r 2 )) 

</R2 



-(Vu, V(Vm,x)) 



1 + r 



2^2 



(W(u),£)ob 



•■J \ (W| 2 + ^(V|Vw| 2 ,x>) + (1+ 2 r2)2 {VV( M ),i) dx 



(6.6) 



(1 + r 2 ) 2 (1 + r 2 ) 3 
2 -l\ „, , 4 



V(u) 



1 + r 2 J y '(1 + r 2 ) 2 
Combining (16. ip . (16. 4p and ( 16. 6 p we obtain 
d 



V{u) dx 
dx. 



dt 



w(r )\u t \ + 



1 + r s 



V(u) 



1 + r 



2^2 



dx 



(6.7) 



l-x 3 



\u t \ 2 + x 3 V(u) dv g 
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□ 



We view Proposition 16.21 as a tool for studying vortex annihilation for fixed e. Note for ex- 
ample, that if all vortices initially reside in the first quadrant so that V(u(x, 0)) is significantly 
larger in the first quadrant than in all others, one would have 

/ XiV(u)dv g > at time t = for i = 1,2, 3. 
Js 2 

Hence, along with the standard energy E e , by Proposition 16. 2\ the three weighted energies 
would also dissipate for some positive interval of time. This could provide the basis for an 
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annihilation result for ( 13. ip analogous to Theorem 15.71 for the ODE flow. Carrying out such 
an argument remains work in progress. 
Finally, we note the following: 

Corollary 6.2. Any stationary solution u to (13. ip . or equivalently, to ( 14. ip on S 2 must satisfy 
the first moment identities 



One can view these conditions as in some sense balance laws for the placement of vortices 
in critical points of Ginzburg-Landau on S 2 . In this light, they can be compared with the 
balance condition derived in Corollary [521 for collections of vortices representing critical points 
of the renormalized energy W (cf. ( 15.11) ) on S 2 . It remains a challenging open problem to 
establish a type of symmetry result for the vortices of critical points on S 2 but (16. 8p at least 
points in this direction. 
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